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ABSTRACT 

Some clusters of galaxies in addition to thermal bremsstrahlung (TB), emit de- 
fNl ■ tectable diffuse radiation from the intercluster medium (ICM) at radio, EUV and hard 

x-ray (HXR) ranges. The radio radiation must be due to synchrotron by relativistic 
tJ" ! electrons, and the inverse Compton (IC) scattering by the cosmic microwave background 

radiation of the same electrons is the most natural source for the HXR and perhaps the 
EUV emissions. However, simple estimates give a weaker magnetic field than that sug- 
gested by Faraday rotation measurements. Consequently, non-thermal bremsstrahlung 
(NTB) and TB have also been suggested as sources of these emissions. 
O We show that NTB cannot be the source of the HXRs (except for a short period) 

and that the difficulty with that the low magnetic field in the IC model is alleviated 
if the effects of observational selection bias, non isotropic pitch angle distribution and 
spectral breaks in the energy distribution of the relativistic electrons are taken into 
account. From these consideration and the strength of the EUV emission, we derive a 
spectrum for the radiating electrons and discuss possible acceleration scenarios for its 



' productions. 



We show that continuous and in situ acceleration in the ICM of the background 
thermal electrons is difficult and requires unreasonably high energy input. Similarly ac- 
celeration of injected relativistic electrons, say by galaxies, seems unreasonable because 
it will give rise to a much flatter spectrum of electrons than required, unless a large 
fraction of energy input is carried away by electrons escaping the ICM, in which case one 
obtains EUV and HXR emissions extending well beyond the boundaries of the diffuse 
radio source. A continuous emission by a cooling spectrum resulting from interaction 
with ICM of electrons accelerated elsewhere also suffers from similar shortcomings. The 
most likely scenario appears to be an episodic injection-acceleration model, whereby one 
obtains a time dependent spectrum that for certain phases of its evolution satisfies all 
the requirements. 
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Subject headings: Galaxies: Clusters: General-Galaxies: Clusters: Individual-X- 
Rays-Magnetic Fields-Acceleration of Particles. 

1. INTRODUCTION 

The most prominent radiation from the intercluster medium (ICM) of clusters of galaxies is 
the thermal bremssstrahlung (TB) or free-free emission in the soft X-ray (2 to 10 keV, SXR) region 
which can reach a luminosity Lsxr ~ 10 45 erg/s and implies gas temperatures of T ~ 10 8 K and 
emission measures of EM ~ 10 68 cm -3 (density n ~ 10~ 3 cm -3 , radius R ~ 1 Mpc). For Coma 
cluster Lsxr — 5 x 10 44 , kT = 8.2 keV. There is, however, a growing evidence for a significant 
nonthermal activity in some clusters. The first of these to be discovered in just a few clusters, 
notably in the Coma cluster (for the most recent observations see Giovannini & Feretti 2000) was 
the diffused (so-called halo) radio emission of luminosity ~ 10 41 erg/s in the frequency range 
30MHz < v < 4GHz whose source must be synchrotron radiation by relativistic electrons. The 
range and distributions of the Lorentz factor 7 of these electrons and their total energy depends 
on the strength, geometry and distribution of the magnetic field B. The field is measured by 
Faraday rotation to be a few microGauss (//G) in some clusters which would require electrons with 
7 > 10 3 . The exact source of these electrons is still a matter of considerable debate. For a review 
see Eilek (1999), Giovannini et al. (1993), Kim et al. (1990), and references there to earlier works. 
More recently, radiation (most likely nonthermal in origin) have been discovered in form of excess 
flux at low and high ends of the thermal radiation in several clusters. The Extreme Ultraviolet 
Explorer (EUVE), with the passband of 69 to 245 eV, has detected excess radiation from Coma 
(Lieu et al. 1996a). The luminosity of this soft excess radiation in the 0.07 to 0.4 keV range is 
Lew ~ 2 x 10 43 erg/s (Lieu et al. 1999). Similarly, BeppoSAX and Rossi X-ray Timing Explorer 
(RXTE) have detected excess hard X-ray (HXR) radiation in the 20 to 80 keV range from Coma 
(Fusco-Femiano et al. 1999, Rephaeli et al. 1999), A2256 (Fusco-Femiano et al.2000), and possibly 
A2199 (Kaastra et al. 1999). The luminosity of Coma in this range is Lhxr ~ 4 x 10 43 erg/s. 
The EGRET instrument on board the Compton Gamma-Ray Observatory (CGRO) set an upper 
limit of Lj-Ray < 10 43 erg/s above 100 MeV (Sreekumar et al. 1996). These observations are 
summarized in Table 1 and assume a Hubble constant of 60 km/(s Mpc). 

Initially the excess EUV radiation was explained in terms of one or more cooler thermal 
components (Lieu et al. 1996a and 1996b, Mittaz et al. 1998) but soon after, a nonthermal 
process, namely the inverse Compton scattering (IC) by the cosmic microwave background (CMB) 
radiation of electrons with energies about ten times smaller than those responsible for the radio 
emission was propsed as the source of the EUV radiation by many authors (Hwang 1997, Enfilin & 
Biermann 1998, Bowyer & Berghdfer 1998, Lieu et al. 1999). The initial interpretation of the HXR 
excess was also based on the IC model (Fusco-Femiano et al. 1999, Sarazin & Lieu 1998). This 
seems to be a natural explanation since electrons with energies very similar to those producing the 
synchrotron emission are required. In fact, long before these discoveries strong upper limit on the 
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nonthermal X-ray emission were set based on the IC model and the radio observations (Schlickeiser 
et al. 1987). 

There is, however, a major difficulty with both of these interpretations (Bowyer & Berghofer 
1998, Fusco-Femiano et al. 1999, Rephaeli et al. 1999). This is due to the simple fact that the 
ratio of the IC to synchrotron luminosities is equal to the ratio of the CMB to magnetic field energy 
densities which for Tqmb = 3K is 15/ (B / fiG) 2 . The observed ratio of HXR to radio luminosities of 
about 4 x 10 2 implies a field strength B < 0.2//G which is much smaller than B values of several fiG 
deduced from Faraday rotation (Eilek 1999) and equipartition of magnetic and relativistic particle 
energies. Comparison of the EUV and radio fluxes can also set a limit on the magnetic field but 
here the limit is somewhat higher (B < 1/xG; Hwang 1997, Enfilin Sz Biermann 1998) and less 
reliable because it is sensitive to the uncertain extrapolation of the electron spectrum over a decade 
(Bowyer & Berghofer 1998). Because of this discrepancy, several workers have proposed nonthermal 
bremsstrahlung as the source of the observed HXRs (Enfilin et al. 1999, Sarazin & Kempner 2000, 
Blasi 2000a). However, this explanation also suffers from a major flaw because it requires a large 
input of energy in the ICM whose consequences have not been detected. This flaw is based on the 
simple fact that bremsstrahlung is an inefficient mechanism. 

In the next section we describe some details of the characteristics of the ICM plasma and the 
constraints they put on the models. In §3 we discuss the emission process and in §4 the related 
aspect of the particle acceleration. A brief summary is presented in §5. 



Table 1 



SOME COMA CLUSTER OBSERVATIONS 



Radiation 



Radio 



EUV 



SXR 



HXR 



Range 
Spectrum 
Index (or Temp.) 

Lum. (erg/s) 
Mechanism 



0.03 - 4 GHz 
Broken Power Law 
~ 1 and ~ 2 
Exponent. Tail 
10 41 
Synchrotron 



0.07 - 0.4 keV 
Uncertain 

~ 3/2 
kT ~ 2keV 
2 x 10 43 
IC (TB) 



2-10 keV 



Exponential 
kT = 8.5tg|keV 
k = 7.51 ±0.18 keV 
5 x 10 44 
TB 



20 - 80 keV 
Power Law 

0.7 to 6 a 
2.35 ± 0.45 6 
4 x 10 43 
IC (NTB) 



a) Fusco-Femiano et al. (1999) b) Rephaeli et al. (1999) 



2. GENERAL DESCRIPTION OF THE PROBLEM 

In order to illustrate the difficulties faced in the above models, in Figure 1 we show the energy 
loss timescales, ti OS s = —E/E\ oss , as a function of particle kinetic energy for all the relevant processes 
in this problem. Here and in what follows unless explicitly expressed all energies and loss rates will 
be in units of rest mass energy of electron, m e c 2 , so that the Lorentz factor 7 = E+l = (1 — /? 2 ) -1 / 2 . 
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We use the following expressions for the loss rates. 



the = (327r/9)r 2 C /3 2 7 \ ph = 2.04 x 1Q~ 20 (3 2 ^(T CMB /3) 4 s" 1 , 

4ych = (4/9)rg C /3 2 7 2 5 2 = 1.32 x 1(T 21 /3 2 7 2 (5 //iG) 2 s" 1 , 
£ Coul = 4vrr 2 cnlnA//3 = 1.20 x 10" 15 (n/10" 3 cm" 3 )//3 s _1 , 



(1) 



(2) 
(3) 
(4) 



£brcm = (16/3)ar 2 cnp 1X (E) = 9.29 x 10- 20 (n/10- 3 cm" 3 )/3 7X ( J E) s 



-i 



where a is the fine structure constant, u p h is the soft photon energy density. The IC losses are 
evaluated assuming the CMB with temperature of 3 K as the source of the soft photons. The 
synchrotron losses are evaluated for an isotropic pitch angle distribution and the Coulomb logarithm 
InA is set to 40, a representive value for the ICM conditions. For the Coulomb and bremsstrahlung 
losses we assume presence of 10% (by number) of fully ionized helium. The bremsstrahlung rate 
also depends on the complex function x(-E) which is equal to one in the nonrelativistic limit and 
equals |[ln(2£ l ) — |] at extreme relativistic energies. (For E S> oT l the slow varying term in the 
square brackets tends to the constant value of about 5.) The bremsstrahlung rate used in Figure 1 
is calculated from the more exact expression given by the Formula 4BN of Koch & Motz (1959). 

Several immediate conclusions can be drawn from the above figure. 

The lifetimes of electrons with energies in the range 200keV < E < 200GeV are longer than the 
free crossing time of the electrons across the cluster (or the 'mean free paths', X\ oss = c/3ti oss < 1 
Gpc, are much larger than the size R ~ 1 Mpc of the cluster). Therefore, these electrons, if 
unhindered, e.g. by chaotic magnetic fields or other scattering agents, will escape the cluster before 
losing most of their energy and while in the cluster they will radiate what is commonly referred 
to as a thin target spectrum. The escaping electrons will radiate most of their energy outside the 
cluster, presumably by IC scattering of the CMB photons. This will disagree with the observations 
and will require a higher rate of energy input than for electrons outside the above energy range 
which lose all their energy before escape and develop a so-called cooling spectrum and giving rise 
to a thick target photon spectrum. Therefore, if all electrons were to lose all their energy in the 
cluster, they must be trapped efficiently so that they traverse a Gpc in the ICM. This can come 
about by a thousand reversals of the magnetic field lines or a million random scattering of the 
electrons. Hence, 

we require the presence of scattering agents (e.g. plasma turbulence) with a scattering mean free 
path A scat ~ 1 pc, or a chaotic field structure with a scale B/AB between 1 kpc to 1 pc. 

In such a case, the electrons with energies away from the peak energy E p of the total loss time 
curve at about 100 MeV diffuse through only a distance of A e gf ~ (Aio^A^t) 1 / 2 <C R before they 
lose most of their energy. This means that for production of a smooth diffuse radiation throughout 
the cluster 

we need in situ acceleration of the electrons throughout the ICM 
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0.01 0.1 1 10 100 1000 10 4 10 5 

Electron Kinetic Energy, E (in mc 2 ) 



Fig. 1. — The energy loss timescales vs energy for the four relevant interactions of electrons for 
typical ICM conditions. The three solid lines, from top to bottom, are for the combined synchrotron 
and inverse Compton, the three radiative processes, and all losses, respectively. The dotted lines 
show the average crossing time T cross ~ R/(c/3) across a region of size R ~ IMpc, the scattering 
time r scat ~ A SC at/(c/3) for a constant scattering mean free path A sca t, and the escape time T esc ~ 
^c 2 ross/ r scat- Note that all loss times are shorter than the Hubble time (heavy solid line for Hubble 
constant of 100 km s _1 Mpc -1 ), and much longer than the electron crossing time except at low 
(E < 200 keV) and very high (E > 200 GeV) energies. 
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and cannot rely on injection of accelerated electron into the ICM from a single source or sources 
separated by a distance 3> A c g. This last condition is required for electrons with E < 200 keV and 
E > 200 GeV even in the absence of scattering agents. (Berezinky, Blasi &; Ptuskin, 1997, address 
the issue of the confinement of the non thermal particles using the scattering mean free path given 
by Schlickeiser et al. , 1987. However, their discussion is applicable only to protons, and other ions, 
because they ignore radiative and other losses. As evident from Figure 1 this cannot be the case 
for electrons.) We will return to these requirements in §4 dealing with the acceleration process. 

On the other hand, because the lifetimes at all energies are much shorter than the age of the 
universe, then unless the observed nonthermal radiations are short lived transient phenomena, 

the acceleration of the nonthermal electrons must be continuous over the lifetime of the clusters, 

Finally, we note that for B < 1/iG, the synchrotron process has little influence on the dynamics 
(acceleration and cooling) of the nonthermal electrons. It acts only as the radio emission process. 



3. EMISSION PROCESSES 

In this section we describe the difficulties faced in some of the proposed radiation mechanisms 
and derive a spectrum for the nonthermal electrons. 



3.1. NonThermal Bremsstrahlung Emission 

EnBlin et al. (1999) were the first to propose this emission (NTB) as the source of the observed 
HXR flux from clusters, whereby electrons of comparable or slightly larger energies produce the 20 
to 80 keV radiation. Sarazin & Kempner (2000) evaluated bremsstrahlung spectra using various ac- 
celerated electron spectra and detailed bremsstrahlung cross sections. Blasi (2000) gives a combined 
description of the stochastic acceleration and bremsstrahlung radiation. However, all these works 
ignore the huge energy problem associated with this model. As is evident from Figure 1, the main 
difficulty of this model is the inefficiency of the bremsstrahlung process compared to the collisional 
losses for E < 1 GeV and relative to IC losses for E > 10 MeV. In particular, for the energy range of 
interest here (20 to ~ 1000 keV) the ratio of the bremsstarhlung to Coulomb loss rates is less than 
10 3 . As shown by Petrosian (1973) the yield of the bremsstrahlung photon is a well defined quantity 
independent of many unknowns of the models. For a non relativistic electron of initial energy Ei n 
that loses all its energy (thick target case) this yield is lb rem = (4/37r)(a/lnA)£'j ra = 7.7 x 10~ 5 Ei n . 
(This yield is larger by a factor of two for electrons losing a small fraction of their energy in 
the source region, thin target case). For a power law distribution of electrons (N(E) oc E~ p , for 
E > Ei n ), the above expression is modified by a factor of order unity: The yield of electrons with 
energies between Ei n and Ef will depend on p. For p ~ 3.5 required by this model this factor is 
1.3 (see eq. [31] in Petrosian 1973, where (5 + 1 = p). This expression is also valid for relativistic 
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energies within a factor of the order of lnEi n as indicated by the slow decline of Tb rcm curve in 
Figure 1 at high energies. 

A yield of Y brcm < 3 x 1(T 6 in the 20 to 80 keV range means that, independent of most details 
of the acceleration or emission model, a large amount of energy (L- m = LuxR/Ybvcm ~ 10 49 ergs/s) 
is fed into the background plasma. If the ICM plasma were to cool only radiatively (free free 
emission), at the very slow rate of L G = 1.45 x 10- 23 ergs/s(T/10 8 K) 1 / 2 ( J EM/10 68 cm °), then such 
an input of energy will increase its temperature at the rate of dT/dt = Li n /3Nk > 10~ 7 K/s. As a 
result the ICM temperature will exceed 10 8 K after a short time of 3 x 10 7 yr and will exceed 10 10 K 
in a Hubble time! This, of course, is not acceptable because it will evaporate the ICM plasma into 
the general intergalactic medium. Either only one part in I0 4or5 of the observed HXR flux is due 
to the NTB process or the NTB emission phase at the observed rate is a short lived phenomena. 
Blasi (2000a) finds that his acceleration model indeed requires a similar (though somewhat smaller) 
rate of input of energy into the turbulence needed for acceleration, and that the duration of NTB 
emission satisfying the observation is around several hundred million years. 

The situation is very similar in what one may call the inverse bremsstrahlung model, whereby 
accelerated protons interacting with the background thermal electrons produce the HXRs. In the 
rest frame of an accelerated proton of energy E p = (m p /m e )E e the process is identical to that of 
bremsstrahlung by accelerated electrons of energy E e . Thus, HXRs of energy 20 to 200 keV can 
be produced by nonthermal protons of energy 40 to 1000 GeV. However, here again, most of the 
proton energy will go into heating the electrons by inelastic Coulomb or Rutherford scattering. In 
addition, the higher energy nonthermal protons may lose some of their energy to tt° production 
which decay into 50 to 100 MeV gamma-rays. 

The presence of thermal SXRs and nonthermal HXRs (also nonthermal synchrotron radio) 
emission in the clusters is very similar to that observed in solar flares. Except in solar (and most 
likely in other stellar) flares the SXR flux is 10 5 to 10 6 times larger than the HXR flux in agreement 
with the above yields. In analogy to flares one may consider acceleration of electrons is taking place 
in high density magnetic loops associated with the disks or halos of, say a thousand galaxies, each 
receiving 10 45 ergs/s. The current observations do not have the spatial resolution to distinguish 
the ICM emission from that of many galaxies. Since the radiative equilibrium temperature T oc 
(L- m / '{Yt, rera EM)) 2 , a lower temperature of about 10 9 K will result for a L m ~ 10 45 ergs/s, density 
n ~ 0.1 and region of size L ~ 30 kpc. (Actually for this size scale conduction losses L con( j ~ 
2 x 10 45 (T/10 8 K) 3 - 5 (L/30kpc) become comparable and exceed the radiative losses for T > 10 8 K so 
that the temperature never exceeds this latter value.) In any case such hot galactic plasmas will 
evaporate into the ICM and may be the source of the hot SXR emitting gas. Only a small fraction 
(< 0.01%) of the energy can go into the ICM plasma. Most of it must be dissipated in the galaxies. 
It is not obvious how the effects of such an energy input which is much larger than that from stellar 
sources (stellar winds, supernovae and other explosions) can be hidden. 

Therefore, we conclude that the main objection to the NTB emission is very robust; it is 
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essentially determined by the values of the fine structure constant and the Coulomb logarithm and 
very difficult to circumvent. This leads to the conclusion that for all three clusters, Coma, A2256 
and A2199, 

the NTB emission from ICM as source of the observed HXRs is not tenable, unless it is a short-lived 
(< 10 8 yr) phenomenon. 

A corollary of this is that one can put a strong constraint on the spectrum and energy density of 
the nonthermal electrons below the peak energy E p ~ lOOMeV where the elastic Coulomb collision 
loses are larger than the radiative losses. As we will show below, the spectral distribution of the 
electrons below this energy must be flatter than E^ 1 ^ 2 or there must exist a sharp cut off below 
several MeV. 

3.2. Inverse Compton Emission 

As can be seen from Figure 1, for typical ICM conditions, the IC emission exceeds bremsstrahlung 
for E > 10 MeV. However, as already pointed out by many of the authors cited in the introduction 
this model also suffers from the inefficiency of the IC radiation relative to the synchrotron radia- 
tion. The relative flux of these two radiations depends on the CMB photon density, which is known 
accurately, and on the value of the magnetic field which is not so well known. From equations (1) 
and (2) the ratio of these fluxes can be obtained to be roughly equal to 

Eic/E sych = 19.8(T CM b/2.8K) 4 (^G/ j B) 2 . (5) 

The most reliable measures of ICM B field come from the Faraday rotation of the background radio 
sources. In the cores of several well studied clusters values of several iiG have been derived (Eilek 
1999). Furthermore, These refer to the net line of sight component so that for a chaotic field the 
actual value could be even larger. This is in apparent contradiction with the value B ~ 0.2^iG 
one deduces from the observed ratio of the HXR to radio fluxes, which is about 500. Actually, the 
observed Faraday rotation in the Coma cluster gives a value of B ~ 0.3/uG for an ordered magnetic 
field and a larger value of B ~ 2//G if the field is chaotic on the scale of several tens of kpc (Kim 
et al. 1990, Ferretti et al. 1995). Estimates based on the assumption of energy equipartition 
between nonthermal electrons, protons and magnetic field give B ~ l^G. The radio properties 
are somewhat different for A2256 but the same kind of discrepancy seem to be present for this 
cluster as well (Fusco-Femiano et al. 2000). However, there is no Faraday rotation data for this 
cluster so the argument against the IC model is not as strong here. The detection of HXR in A2199 
is generally considered as marginal, which combined with the absence of a detectable halo radio 
source or Faraday rotation makes conclusions based on this cluster less reliable. 

However, in comparison with the insurmountable difficulty of the nonthermal bremsstrahlung 
model, there are possible ways to avoid the problems of the IC model. Below we describe several 
effects that alleviate these problems. We will use Coma cluster for quantitative discussion. 
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3.2.1. Selection Effects 

The low value of the B field in clusters with observed HXR emission can be simply an ob- 
servational selection effect. For a given radio flux of the halo source, and independent of any 
equipartition argument, clusters with lower B values will have the stronger IC flux and, therefore, 
will be more readily detected by BeppoSAX and RXTE (or EUVE if electron spectrum extends to 
lower energies). Unfortunately the numbers of known clusters with either (or both) a halo radio 
source and HXR emission are two small to make any reliable quantitative estimates of the effects 
of this selection bias. 

A related and similar effect can arise if the distrubutions of the magnetic field and relativistic 
electrons are inhomogeneous and anticorrelated. In this case the radio and IC emissions will 
come mainly from weak field regions while the Farady rotation is determined by the average field. 
Even in the absence of such an anticorrelation, there are other subtle effects arising from spatial 
inhomogeneities that can give rise to a discrepancy between the magnetic field strengths based on 
the IC emission and the Faraday rotation measure (Goldshmidt & Rephaeli 1993). Explorations of 
spatially inhomogenous models is beyond the scope of this paper. 



3.2.2. Complex Electron Spectra 

The estimate of ratio of IC to synchrotron emission based on equation (5) is for a monoenergetic 
electron. For a spectrum of accelerated electrons this relation is somewhat more complex. However, 
for a power law distribution of accelerated electrons with index p one obtains similar constrains on 
the value of the magnetic field using the observed HXR and radio fluxes. Using the well known 
expression for the spectra of IC and synchrotron emissions it can be shown that the ratio of the 
HXR photon flux (in units of ph/(s cm 2 keV) at photon energy e to the radio flux (in Jy units) at 
frequency v is 

R= J-^) =1MXW Um^keVJy J V^kev) vGhJ U§kJ ) 9 ^ 

(6) 

where a = (p+ l)/2 is the IC photon spectral index and g(p) is a complicated function of index p 
which is equal to 11.0, 41.2, 181 and 755 for p = 2, 3, 4 and 5, respectively (see eqs. [6.36] and [7.29], 
Rybicki & Lightman 1979). In this range of p a good approximation to use is g(p) = e ( 1A2 P~ - 51 ) . 
Using this approximation it can be shown that 

{BJtiG) = (20keV/e)(z^/GHz)f+Texp{2.84(p + r)/(p + 1)}, (7) 

r = 0.71n(i? o6s (e, v)/l.l\ x 1(T 8 ), (8) 

where R b s (e,v) is the observed ratio of the fluxes. Using the observed flux ratios from Coma at 
several values of e and v we find that for p = 3 the field strength is B± = 0.18//G. The required value 
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of -B_|_ increases with p monotonically but slowly. For example, for p = 5 we find B± = 0.8 to 0.3/iG 
depending on the values of e and v. This indicates that magnetic fields of about l^G may be 
possible if the electron spectrum steepens at some energy just below that needed for production of 
HXRs. 

The energy range of electron needed for production of observed HXR (20 to 80 keV), EUV 
(0.07 to 0.4 keV) and radio (0.03 to 3 GHz) are 

0.53 <(E HXR /W 4 )< 1.2, (9) 
0.3 < (E EUV /10 3 ) < 0.75, (10) 
O^G/BjVa < {Erad/w ^ < 4 ( / uG/ J B ± ) 1 /2. (n) 

Note that radio waves with v > 0.35GHz are emitted by electrons above the range needed for the 
other emissions. Thus, a steepening of the accelerated electron spectrum at E = E cr ~ 10 4 will 
reduce the radio flux and allow a higher magnetic field. For example, if the spectral index of the 
electrons changes from 3 to 5 at E cr (as in the Rephaeli 1979 model), then following equations (6) 
and (7), it can be shown that we need B± ~ 0.5/[/G(-E cr /10 4 ) -2 . Even higher magnetic fields will be 
allowed if the electron spectrum cuts off exponentially, as is the case for some acceleration models 
described by Schlickeiser et al. (1987, see also below ). For N(E) = N (E/E cr )-Pexp{E/E cr }, the 
radio flux at high frequencies is reduced approximately by a factor of (p — \ /Z)E p+2 i^[{v / v^) 1 / 2 ], 
where v CI = 0.42GHz(Sj_///G)(-E cr /10 4 ) 2 and E n (x) is the exponential integral function. (This re- 
sult is obtained by approximating the monoenergetic synchrotron spectrum as r](y, E) = A(y /v c ) 1 ^ 
for v < v c = 3E 2 v B /2). Schlickeiser et al. (1987) show that a power law spectrum with an expo- 
nential cut off at u cr = 0.15GHz provides a much better fit than a single or a double power law 
model. With this cut off frequency the application of the above correction factor yields B± ~ 1.7/LtG 
and B± ~ 1.1/iG for p = 3 and 4, respectively. These higher field strengths are in better agreement 
with the Faraday rotation measures quoted above. 

We will return to these consideration in §3.3 and 3.4 and show that these requirements set 
further constraints on the acceleration mechanism. 

3.2.3. Anisotropic Pitch Angle Distribution 

The gyroradius of the nonthermal electrons r g = 2-kc(3j /vb^ ~ 10 11 cm/37(//G/Bj_) is much 
smaller than all other relevant scales in clusters. Therefore, the electrons are attached to the field 
lines and their distribution can be described by a gyrophase averaged distribution N(E, ip) as a 
function of energy and pitch angle ip. In the above discussion we have implicitly assumed an iso- 
topic pitch angle distribution. Anisotopies can modify some of the results quoted above. Note that 
all values of magnetic field are quoted in terms of B±. For an ordered field and isotropic pitch 
angle distribution the synchrotron emissivity is related to the component of the field perpendicular 
to the line of sight. However, as stated above the magnetic field, even though ordered on the scales 
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comparable to and larger than r g , must be chaotic on kpc scale. The emissivity averaged over scales 
larger than one kpc will be isotropic independent of any anisotropics in the monoenergetic emissiv- 
ities and the pitch angle distribution. However the overall intensity will depend on the pitch angle 
distribution. In this case B± = Bsirn^. The synchrotron emissivity at a given frequency is pro- 
portional to _B^ +1 )/ 2 so that the above field values must be corrected by the value of (sin^)^" 1 " 1 )/ 2 
averaged over the pitch angel distribution in the range < ip < tt/2. If the distribution is isotropic 
the average value of this quantity is 2/3 and 8/15 for p = 3 and 7, respectively, so that the actual 
values of magnetic field will be 1.5 to 2 times larger than those quoted above and could be as high 
asB = 3/uG. 

Even higher fields will be required if the pitch angle distribution is anisotropic and is beamed 
along the field lines. For a Gaussian pitch angle distribution of width ipQ < 1 the field strengths 
increase by a factor of ipQ q , where the value of q depends on several factors but is greater than 
one and could be as high as a few. The spectral shape also deviates from the usual power law 
with index a syc h = (p — l)/2 depending on the values of ipo and Yw For further details on this 
see Epstein (1973) and Epstein & Petrosian (1973). Whether the acceleration mechanism will 
accelerate E > 10 4 electrons preferentially along a jumbled field line depends on the conditions in 
the background plasma. We will return to this in §4 wherre we will argue in favor of the isotropic 
distribution. 

In summary, there does not appear to be an insurmountable discrepancy between the field strengths 
required by the IC model for HXRs and the observed values. 



3.3. EUV Emission and Electron Spectral Index 

If the EUV emission is also produced by the IC process the nonthermal electron distribution 
must extend to < 100 MeV. Unfortunately the electron spectral shape in this range is not well 
determined and we must rely on the extrapolation of the spectra from 10 4 MeV range which can 
lead to a large uncertainty. 

As already alluded to in the previous section there has been considerable discussion of the 
radio spectrum and its implication for the electron spectrum. A single power law fit gives a value 
of radio spectral index of about a S ych = 1-5 implying an electron spectral index of p = 4. However, 
as pointed out by Schlickeiser et al. (1987), broken power laws provide a better fit. For example, 
the fit to the Rephaeli (1979) model yields an spectrum with index of about 1 which steepens to 2 
above 0.6 GHz, implying an electron spectral indices of p = 3 and p/i = 5 respectively below and 
above the break energy of E = 1.6 x 10 4 (^G/Sj_) 1 / 2 . Even better fit is obtained for a spectrum 
with an exponential cut off rj(v) oc jy-°- 52 e W°- 15GHz ) ; which means an index of p = 2. However, the 
range of the acceptable low frequency spectral indices is fairly large. For the last model the 90% 
confidence range of a syc h extends from +0.3 to -1.0 implying — oo < p < 3 (see Schlickeiser et al. 
1987). 
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The HXR spectrum is even more uncertain. Fusco-Femiano et al. (1999) give photon index 
0.7 < a < 3.6 which allows 0.4 < p < 6. Rephaeli et al. (1999) give a similar value but a smaller 
range of 1.9 < a < 2.8 so that 2.8 < p < 4.6. The EUV observation when fitted to a power law 
indicate a photon spectral inedex in the range 1.3 to 2.0. If the EUV emission is also due to the 
IC process these values of the photon index indicate a low energy (10 2 < E < 10 3 ) electron index 
in the range 1.6 < p < 3.0. For a summary of these observations see Table 1. 

It therefore appears that a value of p < 3 is consistent with most of the data. A value of p = 3 
implies an IC photon spectrum /(e) oc e~ 2 and equal energy emission per decade. The ratio of the 
observed EUV flux in the 0.07 to 0.4 keV range of 1.5 x 10~ n to the HXR (20 to 80 keV) flux of 
2.2 x 10~ n ergs/(scm~ 3 ) (see Lieu et al. 1999 and Fusco-Famiano et al. 1999) would indicate a 
p ~ 2.9, which is also consistent with the above values. 

In summary, the EUV, HXR and radio data can be fitted by the IC and synchrotron emission in 
a chaotic magnetic field of strength around 1 to 2 [iG, by electrons with the same spectral distribution 
as that needed for the production of the observed radio spectrum via the synchrotron process. 

3.4. Spectrum of Radiating Electrons 

From the above discussions we can constrain the instantaneous spectrum of the radiating 
electrons as follows. We will assume an isotropic pitch angle distribution. 

The radio and HXR observations indicate presence of a power law electron spectrum with an 
index p < 3 and sharp (preferably exponential) cut off at E > E CT ±2 10 4 (/xG/i?_L) 1 / 2 . If the EUV 
emission is also due to IC process the electron spectrum must extend to about 100 MeV with a 
somewhat lower spectral index (p ^2 2.8). At this energy about half of the electron energy is lost 
through Coulomb collisions and about 10% is radiated as bremsstrahlung photons of e < 100 MeV. 
Below E = 200(10 _3 cm~ 3 /n) the collision losses become dominant and for E < 10(10 _3 cm _3 /n) 
bremsstrahlung surpasses all other emissions. Therefore, there is a limit how far this spectrum can 
be extended. It can easily be shown that if the electron spectrum is extended below 20 MeV with 
p = 3 the collisional heating rate of the background thermal plasma will exceed the rate of SXR 
thermal bremsstrahlung emission rate. Since there are other sources of heating of the plasma the 
electron spectrum must cut off rapidly at this or higher energy. However, if the spectral change 
occurs at a higher energy the cut off does not necessarily have to be so severe. For example, a 
spectral break at ^ m i n ~ E p ~ 200 MeV with index pi < 0.5 can be extended to very low energies 
without violating the heating rate threshold. Note that in any case the nonthermal bremsstrahlung 
emission in the X-ray (20 to 80 keV) as well as gamma-ray (> lOMeV) range will be negligible 
compared to the observed HXR flux. It can reach at most about 10% of the total losses at 200 
MeV. Also note that the commonly used spectrum which is a power law in terms of Lorentz factor 
7 (see e.g. Giovannini et al. 2000) has a natural break (flattening) below 0.5MeV. Such a spectrum 
also must flatten much before below 7 ~ 100 to avoid the violation the above mentioned threshold. 
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There has been some discussion (see Blasi 2000b, Bykov, Bloemen & Uvarov 2000, Sarazin 
2000) of a possible constraint imposed by the observed EGRET upper limit of ~ 6 x 10~ 12 ergs/ (cm 2 
s) above 100 MeV (Sreekumar et al. 1996), which is about 0.3 of HXR and 0.6 of EUV fluxes. This 
constraint is not very stringent because the expected bremsstrahlung flux by electrons above this 
energy is about 30 times smaller than their EUV emission or the HXR emission by higher energy 
electrons (see Figure 1). 

In summary, the radiating electron spectrum can be described as 

N(E) = iV n J (E/E CI )-Pe-( E - E ^ if E p < E, for p < 3; 

I (E/E p yP'(E p /E CI )-P if E < E p , for p, < 0.5. 

In the next section we discuss the types of acceleration mechanisms and plasma conditions that 
can give rise to such an spectrum. However, the reader is reminded that all the results in this 
section, including the above equation, are applicable to Coma and only to other clusters with 
similar observational characteristics. 



4. ELECTRON ACCELERATION 

The above spectrum is not necessarily that of the accelerated electrons. It would be in case of 
thin target emission when only a small fraction of energy is lost during the radiation process in the 
ICM, i.e. if n oss > 

T C ross ~ R/c/3. As discussed in connection to Figure 1, this would appear to be 
the case for electrons with 200 keV < E < 200 GeV. However, we face two critical problems if the 
electrons escape the ICM in a time scale shorter than their loss times. The first is that this requires 
an unreasonably high amount of energy for acceleration of the electrons; electrons in the relevant 
energy range radiate less than 1% of their energy in the ICM. The second problem is that escaping 
electrons will continue to produce IC photons outside the ICM and will give rise to EUV and HXR 
emission that extends well beyond the cluster boundary as determined by the TB and radio (halo) 
emissions. This is not what is observed, especially at EUV energies where the source has been 
resolved (see Bowyer & Berghofer 1998). Consequently, all electrons must be trapped, by chaotic 
fields or turbulence, and lose all their energy in the ICM as in a thick target model. It is therefore 
the totality of the acceleration, scattering and loss processes which determine the spectrum of the 
radiating electrons. In this section we discuss some general features of the acceleration process and 
the conditions which can give rise to the spectrum given by equation (12). 



4.1. General Features of Acceleration 

The trapping of the electrons requires that they undergo repeated deflections or scatterings 
such that their effective transport time across the cluster, which we will refer to as the escape time, 
^csc = r cross (i?/A scat ) > Ti oss , where A scat = c/3r scat is the scattering (or deflection) mean free path. 
For this to be true for EUV emitting electrons (E > 200MeV) we need T csc > 3 x 10 9 yr or a 
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-R/Ascat ^ 10 3 ; i.e. we need more than million random scatterings. This implies a mean scattering 
time scale r sca t <3x 10 3 which is more than 10 3 times shorter than the crossing time and much 
shorter than all other relevant times (see Fig. 1). 

A secondary effect of the repeated scatterings is that the pitch angle distribution of the electrons 
will be isotropic. The short mean free path also means that HXR and radio emitting electrons tra- 
verse distances equal to 1/40 and 1 /200th of the cluster radius within their lifetimes. Consequently, 
for a smooth diffuse source the acceleration must be occurring throughout the ICM with inhomo- 
geneity scale smaller than a few kpc, or the resolution of the observation if it is larger. (It should 
be noted that one can impose an ad hoc energy dependence scattering process with mean free path 
Ascat(^) = R 2 1 Aioss(-E') so that the effective range os all electrons is ~ R and T esc (E) = t\ oss (E); 
see below.) We are therefore dealing with essentially a homogeneous and isotropic situation in 
which case the general Fokker-Planck transport equation describing the gyrophase and pitch angle 
averaged spectrum, f(E,t), of the accelerated electrons is simplified to 

% = Wr-^™ - ffiM " - jr^ + Q(E.t). (13) 

Here D{E) and A(E) are the diffusion and systematic acceleration coefficients, Q(E,t) is a source 
term, El is sum of the loss terms given in equations (1) to (4), and T esc (E) is the escape time. 

For stochastic acceleration by turbulence, a second order Fermi acceleration process, 
D(E) = (3 2 D pp and A{E) = ^p-) (§|§ + fE^) describe the diffusive and systematic ac- 
celerations, where D pp is the momentum diffusion coefficient. From these we can define en- 
ergy diffusion and acceleration times r^is = E 2 /D and r ac = E/A. The escape time is related 
to the mean scattering time r scat ~ where D m is the pitch angle diffusion coefficient; 

T esc = T 2 ross D^. For relativistic particles in resonant interaction with Alfven or Whistler waves, 
for example, r ac ~ (/3//3A) 2 7scat> where c(3a is the Alfven velocity (for further details see e.g. Hamil- 
ton & Petrosian 1992). For ICM conditions (5a — 2.3 x 10~ 4 , which means r ac ~ 10 7 r sca t. For an 
efficient acceleration we need an acceleration time which is shorter than both the escape and the 
energy loss times. For the relevant energies this means r ac < 10 s yr, r scat < 10 yr, and T csc > 10 12 
yr. Such a short scattering time may seem unreasonable but is possible. Very roughly this time 
is about (fl e (m e /mp) rl /tiiib) _1 j where Q e ~ 20(B±//j,G) is the electron gyrofrequency, q is the 
spectral index of turbulence energy density, and /turb is the ratio of the total turbulence energy 
density to the energy density of the magnetic field. A Kolmogorov index of 5/3 and /turb ~ 10 -6 
would give a r scat ~ few years but a steeper spectrum will be less efficient and may require unrea- 
sonably large energy density for the turbulence (see also below). Figure 2 shows a comparison of 
these times with the total loss and crossing times from Figure 1. 

The situation is very similar for acceleration by shocks, a first order Fermi process, which 
also requires turbulence for scattering the electrons back and forth across the shock (see e.g. Jones 
1994). In this case we have the additional systematic acceleration term A s h{E) ~ ((3 s h/(3) 2 D^, 
where the shock velocity c(5 s h is of the order of the sound velocity. For ICM conditions (3 s h ~ 3x 10~ 3 
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Fig. 2. — Comparison of the energy dependence of the Total Loss time (heavy solid line, from Fig. 
1) with timescales for scattering (dashed), acceleration (solid) and escape (dot-dashed) of electrons 
due to stochastic (A) and shock (S) acceleration for typical ICM conditions. Two examples are 
given. One with constant acceleration and other timescales, corresponding to an acceleration 
rate A(E) oc E, and a second with variable timescales where A(E) — > a constant at high E, 
corresponding to an exponent q = q' = 1. The dotted line show the average crossing time Across ^ 
R/{c(5) across a region of size R ~ IMpc. The critical energies where the E dependent acceleration 
time is equal to the escape time, and the Coulomb and IC loss times are shown. 
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so that shock acceleration is about hundred times faster than stochastic acceleration. Hence, 
we require a 100 times longer scattering time, which requires correspondingly smaller density of 
turbulence. This and the corresponding escape time are also sketched in Figure 2. 

We should, however, note that in general these time scales are energy dependent. For the 
cases discussed above one expects these time to vary as E 2 ~ q (see Pryadko Sz Petrosian 1997). 
An example of E dependent time scales (with q = 1) are also shown in Figure 2. Note that with 
increasing scattering time the distance diffused by electrons increases and reduces the above men- 
tioned difficulty with the spatial smoothness of the acceleration process. However, for a complete 
removal of this difficulty we need q — 1 and T sca t oc E 2 . 

We now consider several scenarios with opposing and somewhat extreme assumptions. 

4.2. Continuous Acceleration And Steady State Models 

The age of the GHz radio emitting electrons (E ~ 10 4 ) could be as low as 10 8 years (see Figs. 
1 or 2) so that unless the observed nonthermal emission is a short transient event of comparable 
time scale we require a continuous acceleration or injection of nonthermal electrons in the ICM. 
In this case Q{E) is a constant independent of time and if the density and magnetic field change 
slowly, say on a Hubble time scale like the CMB photons, then on this and shorter time scales 
we will be dealing with a time independent or steady state situation with df/dt = 0. Then f(E) 
obtained from equation (13) represents the radiating electrons and must conform to equation (12). 

4-2.1. Acceleration of Thermal Electrons 

The most likely source for the accelerated electrons might appear to be the background hot 
plasma, Q(E) = {yR /2)nE^ /2 \fEe~ E / Eth , where E th = kT/m e c 2 = 0.02. However, this pos- 
sibility suffers from two serious difficulties. The first has to do with the acceleration process. 
Although acceleration by plasma turbulence of low energy (nonrelativistic) electrons is possible 
(Hamilton & Petrosian 1992), the required conditions for it is not the case in the ICM. Presence 
of short wave (or high k vector, k = 2-itub/c ~ 6 x 10~ 10 (B / /zG) cm -1 ) turbulence and a ratio 
of plasma to gyrofrequency of less than one (or Alfven velocity [3a > {m e /m p ) l l 2 ~ 0.023) is re- 
quired (see Pryadko & Petrosian 1997). In the ICM (3a — 3 x 10~ 4 and the value of this ratio is 
100(n/10 _3 cm _3 ) 1 / 2 (/iG/i?) 3> 1. Furthermore, it is not clear how such waves can be excited, and 
even if excited they will be damped quickly because of the high temperature of the ICM (Pryadko 
& Petrosian 1998, 1999). 

The second and more serious difficulty in accelerating the background plasma electrons has 
to do with the high Coulomb losses (already encounterd in §3.1. The acceleration process must 
overcome the heavy losses the electrons will suffer as they are pulled from their low energy state 
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across the energy range 10 keV to several 100 MeVs. In addition, for a reasonable acceleration time 
scale, r ac ~ 10 8 yr, the accelerated electron spectrum will extend into the nonrelativistic region 
with a relatively steep upturn at E < 0.5, where rc ui < r ac (see Hamilton & Petrosian 1992, Park, 
Petrosian & Schwartz 1997). This does not agree with the desired equation (12), requires a high 
level of turbulence (~ 10 48 reg/s, see Blasi 2000a), and will lead to the input of a high amount of 
energy in the ICM as in the NTB model. As discussed in §3.1, this will heat up the ICM plasma 
to above 10 8 K in less than 10 8 yr. 

We, therefore, can conclude that the background thermal electrons cannot be the source for the 
nonthermal electrons, except for a short period of less than 10 s yr. 



4-2.2. Acceleration of Injected Non Thermal Electron 

To overcome both of the above difficulties we require injection of relativistic electrons, presum- 
ably from the cluster galaxies, as the initial source. We first consider the simplest case of a delta 
function injection, Q(E) = Qq5(E — Eq). The acceleration process will distribute these electrons 
above and below Eq and there could be other breaks at critical energies E^ oul , £j£, and E^ c where 
the acceleration time r ac = rcouh^ic and T esc , respectively. Example of these energies are shown 
in Figure 2. For a detailed discussion the reader is referred to Petrosian (1994) and Parks & Pet- 
rosian (1995), and for some examples of complex spectra to Petrosian & Donaghy (1999). Here we 
describe some of the possibilities relevant to the problem at hand. 

Although in certain circumstances the resultant spectrum can be approximated by a power 
law, this is the exception rather than the rule. A power law spectrum over a wide range of energies 
is achieved for simple diffusion coefficients and for negligible loses. For example, for the simple case 
of 

D(E) = VE q ',A{E) = aVE q '-\ and T csc = E S /{9V) (14) 
and for the special case of s = 2 — q' one gets 

\( y E/E ) a - x -V( x2 + t) ) HE>E , 

where x = (a — l + q')/2. Here we follow the notations in Pryadko & Petrosian (1997) and Petrosian 
&: Donaghy (1999), rather than that of Park & Petrosian (1995) and Park et al. (1997) who use q 
for our q' here. In what follows we assume Eq < 200 MeV and concentrate on the spectrum above 

Eq. 

For 9 > a ~ 1, above Eq, N oc E~^ft e \ This is the kind of acceleration model used by 
Schlikeiser et al. (1987). If s ^ 2 — q' the spectrum will deviate from a power law (exponentially 



as in modified Bessel functions and K n ) at the energy E^ c ~ 0\<>~ 2 +<i'\ . This may appear as a 
good explanation for the exponential cut off needed in equation (12) if E^ c ~ 10 4 . However, this 
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would require T csc <C tic below this energy, which as stressed above is ruled out by observations 
and arguments based on energy budget. As shown in Figure 2 for the two acceleration models 
~ T nc/T csc ~ 10~ 2 or 10~ 4 so that E^ c 3> 10 4 . In addition, because acceleration by shocks (if 
these exist in the ICM) is more efficient than by turbulence, the ratio of the systematic acceleration 
rate to the diffusion rate a = (rdifr/r ac ) ~ Ps/Pa ~ 10 2 . In this case, i.e. in the limit — > 0, 
equation (15) reduces to 



so that to obtain the index p = 3 required in equation (12) we need q' = 4. In general q' is 
related to the spectral index q describing the distribution of the wave vector of the turbulence. For 
Alfven waves q' = q so that we require an spectrum of turbulence which is much steeper than the 
commonly assumed value of 5/3 expected for a Kolmogorov spectrum. As described in §4.1 a high 
value of q will require a high level of turbulence especially for the stochastic acceleration model. 

A more reasonable explanation of the required exponential cut off comes from inclusion of the 
losses in equation (13). As mentioned above deviation from a power law is expected at energy E££ ss 
where a specific loss time is equal to the acceleration time. The deviation occurs in the side where 
loss time is shorter. If ti OS s < T ac for E > Ef£ ss (assume to be > Eq), then the spectrum decreases 
sharply (approximately exponentially) above this energy. This situation can arise from IC and 
synchrotron losses, if the acceleration time decrease more slowly than the loss time (oc E~ l ) as in 
the two examples shown in Figure 2. This requires a systematic acceleration rate of A(E) oc E <2 , 
which for scattering by Alfven waves (either in the stochastic or shock acceleration case) implies 
q = q' < 3. As evident from equation (16) this would give rise to an accelerated electron spectral 
index p < 2 which is too small. (For the constant and variable acceleration time scales shown 
in Figure 2, q' = 2 and 1, and p = 1 and 0, respectively.) In the opposite case of q' > 3, the 
losstime ti oss < r ac below E^ ss and the spectrum steepens (becomes softer). This situation clearly 
cannot produce the exponential cut off at high energies and may arise due to Coulomb losses at 
low (perhaps nonrelativistic) energies E < Eq ou1 . For a thorough discussion of all possibilities see 
Park & Petrosian (1995). 

Some of the above discussion is based on analytic solutions which are obtained for simple 
diffusion coefficients. In general these coefficients are more complex (Dung & Petrosian 1994, 
Pryadko & Petrosian 1997) and the resultant electron spectra could have other features such as a 
plateau just before the exponential cut off (Park & Petrosian 1995, Petrosian & Donaghy 1999). 
Testing these more realistic models is beyond the scope of this paper and not warranted by the 
existing observations. 

In summary, we can conclude that an exponential spectral cut off can be produced at E ~ 10 4 
if above this energy either T eS c < Tioss, or T csc > ti oss but r ac < t 1oss . The first possibility is ruled 
out by observations and the second will give rise to a flat electron spectrum with p ~ 1. Although 
the existing radio, HXR and EUV data do not have sufficient spectral resolution to rule out this 
model, a value of p = 1 is barley acceptable. Such a flat spectrum will also exacerbate the problem 




(16) 
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of low required value for B. More importantly, if the EUV radiation is also due to the IC process, 
then the implied photon spectral index of a = (p+ l)/2 = 1 would mean a HXR to EUV flux ratio 
of (80-20)/(0.4-0.07)=200 which is more than two orders of magnitude larger than the observed 
value of less than 2. We note, however, that one can specify a contrived and unphysical energy 
dependence of the acceleration rate which can steepen the spectrum below a GeV to produce more 
EUV photons. We will not discuss such possibilities. 

The above difficulty can not be circumvented even if the injected electron spectrum, instead 
of being narrow as a delta function, is a broad power law; Q(E) oc E~ po for E > E m i n . In this 
case the final spectrum is obtained by the convolution of Q(E' — E) with the above spectra. If we 
use the model of equation (16), for po > q' — 1 this convolution will have no effect above E m \ n and 
the difficulty remains. But in the opposite case, Po < q' — 1, the acceleration process will have a 
negligible effect, and the resultant spectrum will be essentially same as the injected spectrum which 
is now even flatter. Thus, we conclude that 

the steady state acceleration in the ICM of either thermal or non thermal electrons can not produce 
the requisite spectrum for reasonable physical conditions. 

4-2.3. Transport Effects and Cooling Spectra 

Considering the difficulties with the acceleration in the ICM discussed above we now explore the 
possibility that electrons are accelerated somewhere else, presumably in galaxies, and are injected 
into the ICM, where they undergo only scattering and losses. In this case we still need some kind of 
turbulence to scatter and trap the electrons in the ICM, but we assume that these only isotropize 
the electrons and diffuse them spatially but cause neither diffusion in energy nor acceleration. 
As before, the scattering rate determines the escape time in equation (13) where now we set 
D(E) = A(E) = 0. Because we are interested in relativistic electrons, we approximate the loss 
term in equation (13) as 

E L (E)/E P = (1 + (E/E p ) 2 )/t , (17) 

where 

r = (47rrgcnlnA) _1 = 6.3x 10 9 (10 _3 cm _3 /n) yr and E p ~ [(9/8)(nlnA/(up ft +B 2 /8vr))] 1/2 = 235 

(18) 

are approximately the loss time (multiplied by 2) and the energy where the Total Loss curve 
reaches its maximum (see Figs. 1 and 2). Here we have ignored the bremsstrahlung loss and the 
weak dependence on E of Coulomb losses at nonrelativistic energies. Solution of equation (13) 
will then give the effective spectrum of the radiating electrons (referred commonly to as a cooling 
spectrum) that must conform to equation (12). Q(E) represents the average rate of injection of 
accelerated electrons, which in general will have a broad distribution. A reasonable (and convenient) 
form is a power law, Q(E) = Qq(E / E p )~ po with the same low energy constraints as those used 
in connection with equation (12). Because most of the energy of electrons with E < 100 MeV 
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goes into heating the ICM, the rate of injection of energy below this value must be less than the 
SXR thermal luminosity. Therefore, as stated in the previous section the spectrum of the injected 
electrons must drop off sharply below 8 MeV or have an spectral index pi < 0.5. It should also 
be noted that the sources of injection must be sufficiently numerous and have a distribution such 
that they can produce a surface brightness distribution that is as smooth as that observed at radio 
wavelengths. 

It can then be shown that for a finite, and perhaps energy dependent, T esc the steady state 
solution of equation (13) is 



This is a partially cooled spectrum and has a break at x ~ 1 or at energy E cr where T esc = ri oss . For 
x <C 1 or T csc 3> Ti oss one expect a fully cooled spectrum and for the opposite limit, T esc <C ti oss , the 
spectrum is same as the injected spectrum multiplied by T esc . For example, for T csc = T CSC (E / 'E p ) v ~ x 
and for energies above the maximum of the ri oss curve at about 100 MeV, where ti oss oc E^ 1 , a 
power law injected spectrum (for v > and po > 1) gives 



large fraction of the E < E p electrons escape from the ICM, or more accurately from the turbulent 
confining region with a flux of F CSC (E) oc N(E)/T CSC (E). As already pointed out above this is 
in disagreement with the observations. This difficulty is even more severe for a narrow injected 
spectrum, e.g., a delta function. 

For the more reasonable case of T esc 3> ri oss equation (19) reduces to the fully cooled spectrum 
of N(E) = E^ s J*^° Q(E)dE. For a delta function injection at a high energy the spectrum of the 
radiating electrons will vary as E^ s which will be essentially constant up to 200 MeV and then 
decrease with a power law index p = 2. This does not agree with equation (12). For a power law 
injected spectrum N(E) = (Qo/{po — / E p )~ p "ti oss , where ri oss is given by the heavy solid line 
in Figure 2. For po = 2 this will agree roughly with the data but not with the more accurate model 
of equation (12) with a break at E ~ 10 4 . One way to have such a feature is if the injected electrons 
obtain the imprint of the break at their sources. In this case by substitution of equation (12) for 
N in equation (13) with D = A = = we find the necessary injected steady state spectrum 
to be Q(E) oc E 2 'Pe ( ' E / E ^[{2- p)/ E - 1] at high energies and with a similar expression at lower 
energies. This of course is an ad hoc assumption and does not clarify the acceleration mechanism. 

Thus, unless there exists an arbitrary and contrived injected spectrum, we must conclude that 




(19) 
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where E cr = E p (u% sc /tq) 1 / v '. Thus, for po ~ 3 and v ~ and T esc ~ 0.02to we obtain a spectrum 
with a break at E ~ 10 4 , in good agreement with the radio data (Rephaeli 1979 model). However, a 



the steady state injection and cooling model also fails to describe the observations adequately. 
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4.3. Time Dependent Models 

The upshot of the discussion in the previous section is that a steady state acceleration in 
the ICM or modification of a simple accelerated spectra by transport processes in the ICM fail to 
reproduce the general features of the required spectrum. We therefore consider time dependent 
scenarios with time variation shorter than the Hubble time. In this case we consider solutions of 
the time dependent equation (13). We start with the generic model of a prompt one time injection 
of electrons with Q(E,t) = Q(E)5(t — to). More complex temporal behaviors can be obtained 
by the convolution of the injection time profile with the solutions described below. Similar, but 
somewhat different, treatments of the following cases can be found in Sarazin (1999) and Beunetti 
et al. (2000) 



4-3.1. Transport Effects 

We first consider the transport effects in the ICM without any acceleration. In this case the 
time-dependent equation (13), with D = A = 0, T csc independent of time and energy, and El 
constant in time, has the following formal solution: 

f(E,t) = exp{-t/T eK }Q(E'(E,t))E L (E'{E,t))/E L (E), (21) 

where E'(E,t) = T mv (r(E) — t) and r mv is the inverse function of 

roo 

t(E) = / dE/E L (E). (22) 

J E 

Using equation (17) for E L (E), we find t(E)/t = it/2 - tan _1 (.E/.Ep), T inv (x) = cothx and 
E'/E p = (E/E p + tan(£/ro))/(l - (E/E p ) tan(i/ro)), so that for a power law injection, Q(E) = 
Qo{E/E p )~ po , the solution for po > 2 becomes 

M ' ' P1 1 csc ^cos 2 {t/T )[E/Ep + tan(t/T )}Po V ; 

The solid lines in Figure 3 shows the spectral evolution according to this equation for the specified 
parameters and for several times past the injection epoch. At early times the spectrum is a power 
law in the energy range tan(t/ro) < E/E p < l/tan(t/ro), goes to zero at E/E p = l/tan(t/ro), and 
is flat for E/E p < tan(t/Vo). As expected the power law extends to E > 4 x 10 4 needed for radio 
production only for a short period of ro/100 i.e. < 10 8 yr. The power law portion disappears for 
t > 7tto/4 ~ 5 x 10 9 yr and we obtain a degenerate flat spectrum extending to E p . Between this and 
t > tttq/2 ~ 10 10 yr the cut off moves to lower energies and the amplitude drops as tan Po_2 (t/ro). 
In addition, the spectrum decays exponentially on a time scale of T csc so that we need T esc > tq 
which requires presence of turbulence or chaotic fields. Such a turbulence can also accelerate the 
electrons. Therefore, the above spectra are correct if the acceleration time is longer than tq. 
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It is, therefore, clear that either the observable duration of the nonthermal activity in the 
clusters is a rare phenomenon or we need episodic injection of electrons on a timescale of 100 
million years. Whether mergers and resulting shocks, or AGN activities can provide such a source 
is unknown. If this is the case then the rapid cut off at E/E p = 1/ tan(i/ia«o) ma Y mimic the 
exponential form of equation (12), so that with p < 3 this model will be acceptable. 



4-3.2. Acceleration Plus Transport 

A more varied and complex set of spectra can be obtained if we add the effects of diffusion 
and acceleration. Simple analytic solutions for the time dependent case are possible only for special 
cases. Most of the difficulty arises because of the diffusion term which plays a vital role in shaping 
the spectrum for a narrow injection spectrum. For some examples see Park & Petrosian (1996). As 
we have seen for the steady state case the effect of the diffusion is important for a narrow injected 
spectrum. Here we will limit our discussion to a broad initial electron spectrum in which case the 
effects of this term can be ignored. Thus, if we set D(E) = 0, then the solution (21) of equation 
13) can be generalized simply by inclusion of the systematic acceleration term A(E) in El (see eq. 
[14] and [17]) as 

E L (E)/E P = (1 + (E/E p ) 2 - biE/Ep/" 1 )/^, (24) 

where 6 = oT>tqE p = TQ/T ac (E p ) ~ 10 2 or 1 for the shock or stochastic accelerations, respectively, 
and for the parameters described in the previous section. For a general exponent q' one must resort 
to numerical solutions. For the purpose of demonstration of the effects of further ICM acceleration 
we consider the simple case of q' = 2 (corresponding to the constant acceleration timescale of Fig. 
2), which has a solution similar to that shown by equation (23): 

f(F t) - pxn / tlT xn [T + -(E/E p )tMSt/r )/5r- 2 

where 5 2 = 1 — 6 2 /4 and T± = 1 ± 6tan(<5i/To)/(2<5). This solution (valid for b 2 < 4) reduces to 
that in (23) for 6 = 0. For b 2 > 4 we are dealing with an imaginary value for 5 so that tangents and 
cosines become hyperbolic functions with S 2 = 6 2 /4 — 1. For 5 = or b = 2 either form reduces to 

f(E,t) - exp{-i/T csc }Q [E/Ep _ _ i)t/ro]po . (26) 

The dashed line in Figure 3 show the evolution of spectra for this latter case and Figure 4 shows 
the solution according to equation (25) for larger values of b = 5. As expected with acceleration 
one can push the electron spectra to higher levels and extend it to higher energies, but as described 
below this does not significantly alter the above conclusion based on the transport effects alone, 
but improves the situation somewhat. 

For b < 2 the situation is similar to the case 6 = (no acceleration) except that the spectra 
decay more slowly; the degenerate phase of a flat spectrum is reached later and extends to a higher 



-23- 




Elcctron Kinetic Energy, E (in xnc 2 ) 

Fig. 3. — Evolution with time of a power law injected spectrum (top line) subject to Coulomb and 
IC (plus synchrotron) losses as given by equation (23) (solid lines for times t n = 10 n//2 ro,n = —6 
to 0) and with acceleration (b = 2, 6 = 0) obtained from equation (26) (dashed lines for times 
t n = W n / 2 T ,n = -6 to 1). 
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energy. For b > \f(2 + 4/po) ~ 1-83 a local maximum appears during the degenerate phase just 
below the maximum energy (see Figs. 3 and 4). This peak could be very high and narrow. For 
b > 2 the acceleration becomes more and more important and can quickly reverse the decay and 
give rise to a growing spectrum. As before, for early times one gets a power law spectrum at 
tanh(<5t/ro)/(<5T_) < E/E p < T + S / tanh(6t / tq) . The spectrum now can be sustained to a high 
energy for all times: 

E ma jE p = 1 + (b/2)tanh(5t/T )/5/tanh(5t/T )/5 > 5 + 6/2 ~ b. (27) 

Thus, with faster acceleration rate (i.e. b > 50) we can have electron spectra extended above 10 
MeV. However, the period when the spectrum below this energy is a power law is short. The 
degenerate phase is reached quickly when tanh(5t/ro) = y(8/(8 + 2)). For large values of b this 
gives t/r = hi(<5 + 1 + y/{(6 + l) 2 + l))/2<5 -» ln(2£)/(2<5) which is less than 5 x 10 8 yr for b > 50 
implying a short duration for the power law phase. As evident from Figure 4 soon after the electrons 
are reaccelerated to above 10 4 MeV the power law portion disappears. Of course, the situation 
can be improved with a more complex injected spectrum (e.g. a broken power law, see Brunetti et 
al. 2000) or with a time dependent injection and/or acceleration parameters. However, some fine 
tuning may be required to sustain the required spectrum for a period significantly longer than 10 8 
years, which is essentially determined by the temperature of the CMB and the resultant lifetime of 
the E = 10 4 electrons. In any case the additional acceleration in the ICM improves the situation. 




Electron Kinetic Energy, E (in mc 2 ) Electron Kinetic Energy, E (in mc 2 ) 

Fig. 4. — Same as Fig. 3. Left Panel: For b = 5 (solid lines, for t/r = 10~ 3 , 10 -2 , 10" 1 and 1) and 
b = 60 (dashed lines, for t/r c = 10~ 2 - 5 , 10~ 2 , 10~ L5 and lO" 1 ): Right Panel: For b = 2,5,9,17 
and 26 (from lower to upper spectra) and at three different times, t/ro = 10~ L5 , 10 -0 ' 8 and 10~ ' 2 
(solid, dotted and dashed lines, respectively), obtained from equations (25) and (26). 
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5. SUMMARY AND DISCUSSION 

The purpose of this paper has been to investigate the emission mechanisms for the observed 
nonthermal radiation from the ICM of several clusters and to explore possible acceleration scenarios. 
We have used the observations of Coma cluster for our quantitative analyses. The qualitative 
aspects of the results summarized below are quite general, but the specific values of the parameters 
depend on the assumed values of the density, temperature, size, magnetic fiels etc, some of which 
are poorly known and can vary from cluster to cluster. 

For the radiation mechanism, we have come to two important conclusions. 

1) The source of the HXR flux cannot be nonthermal bremsstrahlung emission by semi- 
relativistic electrons because of the extreme inefficiency of this process, unless this is a short-lived 
(< 10 8 yr) phenomenon. 

2) Inverse Compton scattering of relativistic electrons by the CMB photons is a more natural 
process for production of both the HXR and EUV emissions. We have shown that the problems 
with a low value of magnetic field needed for this mechanism (discussed widely in the literature) 
can be alleviated when we include the effects of more realistic (broken power low) spectra and 
anisotropies in the pitch angle distribution of the electrons. Observational selection bias can also 
favor the IC emission at low magnetic fields. 

Combining the requirements of the IC process for HXR and EUV emissions with the require- 
ments of the synchrotron process for the radio emission, we derive a simple spectrum for the 
radiating electrons as described by equation (12). 

Next we investigate the constraints that this spectrum, and other considerations, put on the 
acceleration mechanism. We consider both second order Fermi stochastic acceleration by turbulence 
and first order Fermi acceleration by shocks. We derive parameters for both these mechanisms so 
that they can accelerate electrons to the required energies of E > 10 4 MeV within their life time 
of 10 8 years or shorter. The important conclusions here are the following: 

1) The ICM must contain a high level of turbulence (or other scattering agent) to trap the 
electrons for time periods longer than their loss timescales and much longer than their crossing 
time across the cluster. 

2) Acceleration of the thermal ICM electrons to relativistic energies will be difficult given the 
low value of the Alfven velocity, and more importantly requires input of a large amount of energy 
in the ICM. It will also give rise to an unacceptable spectrum for the IC model. 

3) Steady state acceleration of injected relativistic electrons gives rise to a flatter spectrum 
than desired, or to a HXR and EUV source that extends well beyond the boundaries of the radio 
source. 

4) Steady state cooling of a power law injected spectrum also suffers from the same shortcoming 
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or must involve ad hoc assumptions. 

5) Time dependent models fair much better. A power law injected spectrum, under the in- 
fluence of transport effects alone, can evolve into one with a high energy cut off at E cr ~ fO 4 
(as required by the observations) after a time equal to the energy loss time at this energy, which 
is about 10 8 yr. For later times (t > To ~ 6 x 10 9 yr) the cut off moves to lower energies and 
the spectrum becomes flat below it. If one adds an acceleration agent then the spectrum can be 
maintained above the desired energy for a longer period. This requires an acceleration time scale 
that is shorter than 10 8 yr. But, at such high acceleration rates, the spectrum below this cut off 
becomes flat in a shorter period of time, t ~ rolnfe/fe, where b ~ ro/r ac . This can yield an acceptable 
spectrum for a period of about 5 x 10 s yr. 

The above results mean that either the nonthermal emissions from the ICM are short lived and 
rare events or there is episodic injection of power law spectrum of relativistic electrons on a time 
scale of about I0 8 yr. This, however, still leaves the initial mechanism of the electron acceleration 
unresolved. A likely scenario is that episodic mergers of sub clusters or encounters between galaxies 
can give rise to shocks and turbulence. The initial acceleration can take place in these shocks. The 
spectrum of radiating electrons is a result of transport and further acceleration (by turbulence) 
in the ICM. In such a situation, however, one would expect a different spatial distribution for the 
EUV emission than for HXR and radio emissions. The latter emitted by higher energy, shorter lived 
electrons will be more concentrated around the initial source. Similarly, a radial variation of the 
magnetic field could result in a more (or less) centrally concentrated synchrotron (radio) emission 
compared to the IC (HXR and perhaps EUV) emision. Density variations will affect mainly the 
bremsstrahlung emission relative to the other radiative processes but not the arguments based on 
the bremsstrahlung yield. Temperature variations can effect the spectrum of the turbulence. 

An exact evaluation of the relative spatial distributions at different energy bands will require 
solution of an inhomogeneons Fokker-Planck equation which in turn requires knowledge of energy 
and spatial dependences of scattering and escape processes, as well as spatial variation of density 
and magnetic field. This is beyond the scope of the present paper and not warranted by the existing 
observations. Higher spatial resolution observation will be helpful here. 

Alternative sites of the initial acceleration may be in galaxies, in which case the homogeneous 
model will be a good approximation. However, in this case, in addition to electrons one would 
expect a larger energy input in form of protons. It is likely that protons may be the source of the 
turbulence which is essential for any viable model of nonthermal emission from ICM. 

This work was started while I was a visitor at the Institute For Advanced Studies and Bochum 
University. I would like to thank the support of both Institutions, their staff and acknowledge 
stimulating discussions with Drs. Bahcall, Kumar, Schlickeiser. I would also like to thank Drs. 
Blasi, Bykov, Eilek, Enfilin and Rephaeli for helpful discussions and several suggestions which led 
to the improvement of this paper. 
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